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ABSTRACT 

Using data from the HST Medium Deep Survey, a long-term Key Project, together 
with generically similar archived data, we have discovered evidence for weak gravita- 
tional 'shear' of background field galaxies (I = 22 - 26) in the vicinity of isolated, 
foreground galaxies (1=15- 22), especially those of early type. The statistical lens- 
ing is demonstrated by the slight preferential orientation of the major axes of the 
background galaxies in directions at right angles to the lines joining them to the fore- 
ground ones. The detected shear is at the level of 8(f) = 2.5° ± 1.2° averaged over 
the background galaxies at about ~ 10 half-light radii from 400 foreground ellipticals, 
corresponding to a polarization of about 0.04. A positive but less significant result is 
also reported for foreground spirals. This result excludes the de Vaucouleurs model 
for the mass distribution of ellipticals, but is consistent with elliptical galaxies having 
total M/L ~ 100 and extensive dark halos. 
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1 INTRODUCTION 

The first attempt to detect the weak gravitational shear 
of background galaxies by isolated foreground field galax- 
ies was made by Tyson et al (1984) who catalogued field 
galaxies between mj — 19 - 21.5 (where the J magnitude 
is approximately equal to Johnson B) and searched for po- 
larisation in nearby field galaxies with magnitudes between 
m.j = 22.5 - 23.5. These authors inferred a 2a upper limit 
on the average mass of the brighter foreground galaxies at 
1.7 x10 11 Mq, with a corresponding limit on the rotational 
velocity of 170 km s _1 for a mass envelope extending to be- 
yond 50 kpc (this limit was subsequently revised by Kovner 
and Milgrom (1987) to about 330 km s _1 ). Tyson et al used 
photographic plates in moderate seeing (2" — 3"), and the re- 
sults from these early studies were negative or inconclusive. 
Using CCD data with good atmospheric seeing, however, 
Brainerd et al (1996) have claimed that polarisation has in- 
deed been detected, for background field galaxies between 
23 < r < 24 mag and within about 35 arcs of foreground 
galaxies between 19 < r < 23 mags. 

The Wide Field and Planetary Camera 2 (WFPC2 - 
Trauger et al 1994) on the Hubble Space Telescope produces 
(nearly) diffraction-limited images which are ideal for these 
studies, since ellipticities and position angles can be mea- 
sured routinely to I = 24 in summed exposures of about 30 
mins. on each of 2 - 3 HST earth-orbits. The field of view is 
small (a total WFC area of 5.0 sq. arc mins.) so that data are 
needed from large numbers of fields in order to accumulate 
a large enough number of galaxies to form useful statistics. 
The HST Medium Deep Survey using WFPC2 commenced 
in January 1994 and has since accumulated a sufficient num- 
ber of fields that we have been able to initiate a search for 



these effects. Furthermore, a primary WFPC2 survey con- 
sisting of 28 contiguous fields was performed by Groth et 
al (1994) in April 1994 and these data have been obtained 
from the HST archive. 



2 OBSERVATIONS 

The MDS uses data taken with WFPC2 operating in parallel 
mode, using the F606W filter for the V-band and F814W for 
the Fband (Griffiths et al 1994a, 1994b). We have excluded 
fields dominated by stellar images. The MDS sample is uni- 
formly distributed over the sky with |6| > 20deg, with the 
detection limit of each field ranging over I ~ 23 ~ 25 and 
V ~ 24 ~ 25.5 (see Neuschaefer et al. 1996). For the pur- 
poses of the present study, the WFPC2 MDS dataset has 
been combined with the strip survey of Groth et al. This 42 
arc min long strip consists of 28 contiguous WFPC2 fields 
which extend along the direction of galactic latitude, with 
the survey center at b = +60°. 25 and I = 96°. 35 The total 
area is about 120 square arc minutes when we exclude the 
0.6 sq arc min overlap between adjacent fields. The twenty- 
eight fields were observed for 2800s (4x700s) in F606W and 
4400s (4scll00s) in F814W. 

The calibrated data from the MDS and Groth strip were 
run through the automated object finder and classifier de- 
veloped for the MDS project (Ratnatunga et al 1994). We 
have corrected the geometric distortion near the boundary 
of CCDs using the solution of Holtzman et al (1995) in or- 
der to avoid spurious detection of signal (See Section 5). For 
each object detected, the observed image is fitted with sim- 
ple model profiles (point source, r 1 ^ 4 profile or 'bulge', and 
exponential profile or 'disk') using a 2-dimensional technique 
of maximum likelihood estimation or MLE (Ratnatunga et 
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al 1994). MLE object classification is then based on the lu- 
minosity profile. Structural parameters such as axis ratio, 
size (r M , angular half light radius) and total magnitude re- 
sult from the best model fit, which also gives the position 
angle of the major axis of each galaxy. Near the complete- 
ness level of the survey, the noise from the MLE fits is three 
times less than the corresponding noise resulting from the 
use of image moments measured using standard data anal- 
ysis packages such as FOCAS (Jarvis and Tyson 1981). 

We detected 475+/-25 objects in each of the Groth strip 
fields, and a similar number in MDS fields of comparable 
depth. After removal of stars and duplicate (confused) ob- 
jects we thus produced a catalog of about 28,000 objects, 
of which about 10,000 galaxies are brighter than the com- 
pleteness magnitude of 1=25. Full details of the catalog and 
statistical analyses of the results therefrom are described in 
other publications by the MDS group. 

The foreground 'lens' galaxy candidates are objects 
classified as bulge, disk or (unclassified) galaxies with 15 < 
I < 22 in the object catalogs (1600 total, of which about 400 
are E/S0 galaxies). The high resolution HST data have made 
it possible to separate these potential lensing populations, in 
a way which cannot be done using ground-based data. The 
background 'source' galaxy candidates are objects classified 
as bulge, disk, or (unclassified) galaxies with 22 < I < 26 
(14,000 total). 

3 GALAXY MASS DISTRIBUTION MODELS 

The galaxy mass density profile is modeled as a softened 
isothermal sphere (see Krauss & White 1992; Kochanek 
1996), the shape of which is determined by the core radius 
(r c ), velocity dispersion (it), and truncation radius rt, be- 
yond which the mass density becomes zero. For E/S0 galax- 
ies, the core radius is taken to be r c = 0.05rhl (Kochanek 
1996), and for spirals r c — 0.5r w (Binney & Tremaine 1987), 
where the half light radius r^i is a function of the lumi- 
nosity following the relation described in Im et al. (1995a). 
The velocity dispersion scales as a function of the luminos- 
ity, following the Faber- Jackson relation (for E/SOs) or the 
Tully-Fisher relation (for Spirals; a = V c /^/(2), here V c is 
circular velocity), and we adopt parameters for the relation 
those used by Kochanek (1996) and Brainerd et al. (1996). 
The truncation radius r t is treated as a free parameter which 
is assumed to scale with r^i- The truncation radius deter- 
mines the extent of the dark halo in galaxies, and is thus an 
important parameter for the measurement of the total mass 
to light ratio. 

For elliptical galaxies, we also modeled the mass pro- 
file using the de Vaucouleurs profile, with a constant mass 
to light (M/L) ratio, since the light profile of ellipticals fol- 
low this model (de Vaucouleurs 1948). For the M/L ratio of 
this model, we assume the mass to near infrared light ratio 
of (M/Li)/(Mq/Lq j) — 20h, which roughly corresponds 
to an upper limit of the M/L ratio at the core of elliptical 
galaxies (e.g, de Zeeuw & Franx 1991). For the redshift dis- 
tribution of galaxies, we use the functional form described 
by Efstathiou (1995), controlled by one input parameter, 
the median redshift z me d- The redshift distribution of the 
foreground galaxies is rather well defined by spectroscopic 
surveys (Lilly et al 1995) and by our own data. For our ellip- 
tical galaxies, we obtained z me d = 0.6, estimating redshifts 



photometrically (Im et al. 1996). Lilly et al. (1995) found 
Zmed — 0.5 for galaxies with / < 22, and we will adopt this 
value for spiral galaxies. The value of z me d for the back- 
ground galaxies (22 < I < 26) is unknown but very likely 
Zmed < 3 (Guhathakurta et al. 1990). We will set z me d of 
the background galaxies equal to 2, and explore the effect 
of smaller values. For the observed axis ratio distribution 
of background galaxies, we use the axis ratio distribution of 
disk galaxies in Im et al. (1995b). 

4 DETECTION OF WEAK SHEAR 

The gravitational shear of background galaxy images 
causes them to be preferentially oriented in the direction 
perpendicular to the radius vector to the foreground galaxy. 
The "relative position angle" <f> is the position angle of the 
major axis of the background galaxy, measured counter- 
clockwise from the radius vector to the foreground galaxy. 
When we fitted galaxy models to the images, we implicitly 
assumed that galaxies are symmetric with respect to a 180° 
rotation, and <j> is therefore defined between 0° and 180°. All 
values of <p are in principle equivalent, on the assumption 
that the foreground and background galaxies are physically 
unrelated. Gravitational shear can in principle be detected 
as a deviation of the distribution of 4> from uniformity. This 
deviation depends on the strength of the gravitational shear- 
ing effect, measured by the "image polarization" p, which is 
the induced ellipticity in an intrinsically circular image. In 
the limit of weak shear that we consider, the polarization 
p equals the difference between the tangential and the ra- 
dial stretch due to the lensing. In practice, the gravitational 
shear will lead to a small bias towards values of <j) around 
90°. 

Different measures have been suggested for the devi- 
ation of (f> from a uniform distribution. The most obvious 
method is to consider the "reduced position angle" <j} re d in 
the range (0°, 90°), obtained by taking either <j) or 180° — 4> 
as required; this is equivalent to merging together clockwise 
and counter-clockwise deviations. The average value of the 
reduced position angle is a measure of the non- uniformity of 
the distribution induced by the gravitational shear. A sec- 
ond possibility is to express the distribution as a Fourier 
series; for the range (0°, 180°), only cos(m0) terms are re- 
quired. The effect of the gravitational shear will manifest 
itself primarily as a negative cos(20) term. 

A third possibility is to estimate directly the polar- 
ization as p ?s — (1 — r 2 )/(l + r 2 ) cos (20), where r is the 
measured image axis ratio; this is closely related to the pa- 
rameter x used by Brainerd et al (1995). Use of this third 
indicator is made difficult by the fact that only 30% of the 
(usually faint) source galaxies have reliable measurements 
of the axis ratio. Since the results are comparable to those 
of the Fourier decomposition, this method is not considered 
further here. 

4.1 Results 

Figure 1 shows the average value of the reduced relative po- 
sition angle, as a function of 8/rhi- The average is larger 
than 45° because gravitational shear makes the tangential 
alignment (0 re d ~ 90°) more likely than the radial align- 
ment (4>red ~ 0°). The results of the model calculations are 
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Fig. 1(a). The mean reduced position angle < (j> re d > of source 
galaxies binned as a function of radial distance from the fore- 
ground ellipticals, where the radial distances are normalized by 
the rhi of the ellipticals 
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than 10 half-light radii (i; 30 kpc for a L» galaxy). The ef- 
fect of changing z me d for the redshift distribution is such 
that a stronger signal is obtained using a high z med . But 
reducing z me d does not reconcile the data with the de Vau- 
couleur model. Since early-type galaxies are known to cluster 
more than spirals (Neuschaefer et al. 1996 have shown this 
for the same data set as that used here), several tests were 
made to ensure that the lensing signal was not produced 
by clusters rather than the ellipticals themselves. The fields 
each contain several E/SOs and the tests involved shifting 
the assumed centers of the potentials to points which might 
likely be close to the possible 'cluster' centers. All such tests 
produced null results for the lensing signal. 

For the rest of the foreground galaxies (mostly spirals) , 
the observed signal is considerably weaker (Fig.l-b). This is 
an expected trend because the late-type galaxies have less 
mass and shallower profile near the center than ellipticals. 
But again we see that the observed signal is well matched 
by the isothermal sphere model with rt = 100 rhi- No signal 
is observed in a control sample with stars as "foreground" 
objects (Fig.l-c; See also section 5). 

It is interesting to note that the observed signal shows 
r t >> lOrw. If r t ~ 100 r hi, then the M/L of the lensing 
galaxies is about a few hundred. This value is comparable to 
the M/L of clusters of galaxies, suggesting that 1) the mass 
of clusters can be explained by the sum of the masses of 
cluster members; and that 2) O mci ti er ~ 0.1. A full compari- 
son between the models and the data, with a full discussion 
of the image polarisation, will be presented by Kneib et al 
(in preparation). 



Fig. 1(b). The mean reduced position angle < re£ j > of source 
galaxies binned as a function of radial distance from the fore- 
ground spirals, where the radial distances are normalized by the 
rhl of the spirals 
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Fig. 1(c). The mean reduced position angle < re£ j > of "source" 
galaxies binned as a function of radial distance from the fore- 
ground stars, where the stars have been taken as a control sample 
(8 hi = 0'.'5 is used for this figure) 

compared to the detected weak lensing signal in Fig. 1. The 
observed signal for ellipticals stays strong at least out to 
about 50 half-light radii, and is well matched by a model 
with r t = 100 r h i ■ The de Vaucouleur model and the isother- 
mal sphere model with r t = 10 rhi both fail to match the 
observed signal, suggesting that the M/L of ellipticals is not 
constant and that the mean truncation radius is much larger 



4.2 Fourier-series decomposition 

The results of a Fourier-series decomposition of the distribu- 
tion of relative position angle (f> for the full sample are shown 
in Fig. 2 as a function of /rhi, where 8 is the apparent an- 
gular separation between source and lens centers and rhi is 
the half-light radius of the lens. For this analysis, subsample 
of galaxies more than 150 pixels from the edge of CCDs is 
used in order to avoid false signal due to the geometric dis- 
tortion at the boundary of the CCD (See section 5). Fourier 
terms are normalized as a fraction of the m = term; for- 
mally, the m th term is given by A m = 2(cos(ra0i)), where 
the index i spans all the sources considered. Only the m = 2 
term shows a significant deviation from zero, as predicted 
by the gravitational lensing effect. The relationship between 
polarization and A2 depends on the axis ratio distribution of 
source galaxies; for an average ellipticity of 0.56, as observed 
in our data for galaxies in the range 22 < I < 26 mag, we 
obtain p w 0.26^2- 

Different subdivisions of the sample, by properties of the 
lens and of the background sources, show the polarization 
signal to different degrees. In Fig. 2 we show the amplitude 
of the A2 component separately for ellipticals and spirals, 
as well as for all lenses. The signal is clearly stronger for 
elliptical galaxies than for spirals. (In fact, the signal for 
faint lenses and for spirals is not statistically significant.) 
This is in line with the predictions of the isothermal sphere 
models, shown as thin lines in each plot. 
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Fig. 2. Weak lensing signal as measured by a Fourier-series de- 
composition method. The model curve corresponds to the solid 
line in fig. la 



when using stars as centers (Figure lc). Finally, the geo- 
metric distortion has been modeled and corrected using the 
solution of Holtzman et al (1995); this allows lenses and 
sources from different camera heads to be considered as an 
ensemble, and has been used in the average angle calcula- 
tion (Fig. 1). However, we find that the geometric correction 
does not influence the signal significantly. 

6 SUMMARY AND CONCLUSIONS 

We have discovered evidence for weak galaxy-galaxy 
shear in data from the HST Medium Deep Survey and 
archived HST survey data. The weak lensing signal ex- 
tends out to 50 half light radii from the foreground galax- 
ies, and results in a mean position angle of about 47.5 de- 
grees at around 10 half-light radii for ellipticals, and 46 de- 
grees for spirals (compared with a null of 45 degs.). This 
result is consistent with a simple lensing model for the fore- 
ground galaxies which employs a softened singular isother- 
mal sphere with r t = lOOr^, or at least r t » lOr^;, and 
a* = 220km/ sec for the central velocity dispersion of L* 
ellipticals and V c * = 220km/ sec for the circular velocity of 
L* spirals. However, a constant M/L de Vaucouleurs model 
for the mass profiles fails to fit the observed lensing signals 
for ellipticals and is excluded. 



5 TESTS FOR SYSTEMATIC ERRORS 

The effect of gravitational shear is a subtle one, and there- 
fore special care must be exercised to identify and prevent a 
number of small problems that could give a false signal, or 
mask and/or weaken an existing signal. One such problem 
is the possibility of faint pieces of the lens galaxy itself, such 
as pieces of spiral arms at large distance from the center, 
masquerading as background galaxies . This is potentially 
a problem at separations of 1-2 r M (it may be the origin 
of some of the signal seen at small separation in the whole 
sample and the spirals in Fig. 2), but it is extremely un- 
likely to be significant at 10 r^i, where a strong signal can 
still be seen. Also, it would affect preferentially the spiral 
galaxies, and the observed signal at large radial distances 
is much stronger for ellipticals. The observed signal could 
potentially be affected by other galaxies on the lines of sight 
to the background sources. We have examined those galax- 
ies within separations of lSr^i of the ellipticals, where we 
find that there are fainter ellipticals (by up to 1 mag.) in 
less than 10% of cases; we conclude that the effect on the 
observed signal is negligible. 

Another potential problem is the geometric distortion 
of the WFPC2 field of view. For systems near the edges of 
the field, this induces a small artificial rotation of the galaxy 
images that could result in a spurious alignment. We have 
tested this in several ways. First, a subsample using only sys- 
tems more than 150 pixels from the edges shows the same 
qualitative distortion as the full sample. Second, a test has 
been carried out in which "foreground" and "background" 
galaxies close in pixel position, but in different fields, have 
been associated; any geometric distortion effect would have 
been visible in this sample, but a null signal was detected 
(control sample in Fig. 2). A null signal was also detected 
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